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Abstract. - We demonstrate, that the main universal features of the low temperature experi- 
mental H — T phase diagram of CeCoIns and other heavy- fermion metals can be well explained 
using Landau paradigm of quasiparticles. The main point of our theory is that above quasi- 
particles form so-called fermion-condensate state, achieved by a fermion condensation quantum 
phase transition (FCQPT). When a heavy fermion liquid undergoes FCQPT, the fluctuations 
accompanying above quantum critical point are strongly suppressed and cannot destroy the 
quasiparticles. The comparison of our theoretical results with experimental data on CeCoIns 
have shown that the electronic system of above substance provides a unique opportunity to 
study the relationship between quasiparticles properties and non-Fermi liquid behavior. 



Although much theoretical efforts have been devoted to understand the non-Fermi liquid 
(NFL) behavior of heavy fermion (HF) metals using the concept of quantum critical points, 
the problem is still far from its complete understanding since the experimental systems display 
serious discrepancies with the theoretical predictions [1] . Common belief is that a quantum 
critical point (QCP) is the point where a second order phase transition occurs at temperature 
T — > 0, and where both thermal and quantum fluctuations are present destroying quasiparti- 
cles and generating a new regime around the point of instability between two stable phases [2] . 
Recent experimental studies of the CeCoIns HF metal provide valuable information about the 
NFL behavior near possible QCP due to its excellent tunability by a pressure P and/or a 
magnetic field H [3-6]. The experimental studies have shown that besides a complicated 
H — T phase diagram, the normal and superconducting properties around the QCP exhibit 
various anomalies. One of them is power (in both T and H) variation of the resistivity and 
heat transport [3-8], inherent to both NFL and Landau Fermi liquid (LFL) regimes. The 
other one is a continuous magnetic field evolution of a superconductive phase transition from 
the second order to the first one [9,10]. Above anomalous power laws can be hardly accounted 
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for within scenarios based on the QCP occurrence with quantum and thermal fluctuations. 
For example, the divergence of the normal-state thermal expansion coefficient, a/T is stronger 
than that in the 3D itinerant spin-density-wave (SDW) theory, but weaker than that in the 
2D SDW picture [11]. This brings the question of whether the fluctuations are responsible 
for the observed behavior, and if they are not, what kind of physics determines the above 
anomalies? Fortunately, the direct observations of quasiparticles in CeCoIns, Celrlns, and 
Sr 3 Ru207 have been reported recently [4,12,13]. However, if the quasiparticles do exist, why 
they are not suppressed by the fluctuations? Moreover, these recent facts contradict strongly 
to the theoretical investigations, where quantum phase transitions responsible for the NFL 
behavior considered to be of a second kind so that the quasiparticles arc inhibited near these 
phase transitions [1,2]. 

In this letter we show that these problems can be resolved within Landau quasiparticle 
picture providing that quasiparticles form the so-called fermion-condensate (FC) state [14] 
emerging behind the fermion condensation quantum phase transition (FCQPT) [15]. We show 
that near FCQPT the fluctuations are strongly suppressed while quasiparticles are "protected" 
from above fluctuations by the first order phase transition. We analyze the experimental H — T 
phase diagram of CeCoIns and show that its main universal features can be well understood 
within the theory based on FCQPT. We demonstrate that the electronic system of CeCoIn 5 
can be shifted from the ordered to disordered side of FCQPT by a magnetic field; therefore 
giving a unique possibility to study the relationship between quasiparticles and NFL behavior. 

To study the low temperature universal features of HF metals, we use the notion of HF 
liquid in order to avoid the complications associated with the crystalline anisotropy of solids. 
This is possible since we consider the (universal) behavior related to the power-law diver- 
gences of observable variables like the effective mass, thermal expansion coefficient etc. These 
divergences are determined by small (as compared to those from unit cell of a correspond- 
ing reciprocal lattice) momenta transfer so that the contribution from larger momenta can 
be safely ignored. Let us consider HF liquid characterized by the effective mass M* . Upon 
applying the well-known equation, we can relate M* to the bare electron mass M [16, 17] 
M* = M/(l — NoF 1 (pf 7 pf)/3). Here No is the density of states of a free electron gas, pp is 
Fermi momentum, and F 1 (pf,Pf) is the p-wave component of Landau interaction amplitude. 
Since LFL theory implies the number density in the form x = p f /3tt 2 , we can rewrite the 
amplitude as F 1 (pf,Pf) = F 1 (x). When at some x — xpe, F 1 (x) achieves some critical 
value, the denominator tends to zero so that the effective mass diverges at T = 0. Beyond 
the critical point xfc the denominator becomes negative making the effective mass negative. 
To avoid physically meaningless states with M* < 0, the system undergoes FCQPT with 
FC formation in the critical point x = xfc- Therefore, behind the critical point xfc the 
quasiparticle spectrum is flat, e(p) = /i, in some region pi < p < pf of momenta, while the 
corresponding occupation number no(p) varies continuously from 1 to 0, < no(p) < 1 [14]. 
Here p, is a chemical potential. To investigate the FC state at T = 0, we apply weak BCS-like 
interaction with the coupling constant g and see what happens with both the superconducting 
gap A and the superconducting order parameter «(p) as g — > 0. Let us write the usual pair 
of equations for the Green's functions F + (p,u) and G(p, ui) (see e.g. ref. [16]) 




; G 



u 2 (p) v 2 (p) 



(1) 



(cj - E(p) + i0){u> + E(p) 
where E\p) = £ 2 (p) + A 2 , £(p) = e (p) 



iO) 



Lu-E(p)+i0 w + £(p)-i0' 



H, and the superconducting gap, 




(2) 
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Here v 2 (p) = (1 — £(p)/i£(p))/2, v 2 (p) + u 2 (p) = 1, and simple transformations give 
with k(p) — u(p)v(p). Next we observe from eqs. (j2J and 10) that 

(4) 

It follows from eqs. @, © an d ijU that when g — > the superconducting gap A — > 0, 
while £ = and the dispersion s(p) becomes flat, providing that k(p) is finite in some region 
Pi < P < Vfi making 5 finite. Thus, in the state with FC A can vanish while parameters k(p) 
and S are finite. Taking into account eqs. (J2J) and J3J we represent eqs. JU as follows 

F+= «(p) «(p) _ c _ u 2 (p) v 2 (p) 

uj-E(p)+iO uj + E(p)-iO' uj-E(p)+iO oj + E(p) - iO' [) 

It is directly seen from eqs. (J5J that in the FC state at g — > 0, the equations for functions 
F + (p,(i;) and G(p, ui) take the following form in the region where fc(p) 7^ 



F+(p,w) = -«(p) 



1 1 



G(P,^^ + ^B) (6) 



oj + iO oj — iO 

Here, the factors v 2 (p), u 2 (p) = 1 — v 2 (p) are determined by the condition e(p) = /j, when 
Pi < P < ?>/• Upon integrating G(p, w) over oj we obtain that v 2 (p) — n(p), where n(p) is 
the quasiparticles distribution function. Taking into account the well-known Landau equation 
SE[n(p)]/5n(p) = e(p), we observe that the equation determining n(p) takes the form [14] 

SE[n(p)] 

where E[n(p)] is Landau functional [16]. Equation Q describes the state with FC character- 
ized by the superconducting order parameter «o(p) = ■\/ rt o(p)(l n o(p) where the functions 
n (p) are solutions of eq. (J7J). It is instructive to construct F + (p,oj) and G(p,oj) when g is 
finite but small so that the functions v 2 (p) and k(p) can be approximated by the solutions of 
eq. 0. In that case, S, A and -E'(p) are given by eqs. (0J, J2J) and JSJl respectively. Inserting 
these into eqs. © we obtain functions F + (p,oj) and G(p, w). It is seen from eq. © that A 
is a linear function of the coupling constant g. Since the transition temperature T c ~ A tends 
to zero along with g — > 0, the order parameter «(p) of the FC state vanishes at any finite tem- 
perature so that at T > the quasiparticle occupation number is given by the Fermi-Dirac 
function which we represent in the form e(p,T) ~ fi(T) — Tln{(l — no(p, T))/(tiq(p, T))}. 
Observing that at T — > the distribution function satisfies the inequality < rio(p) < 1 at 
Pi < P < P/, we conclude that the logarithm is finite (therefore Tln(...) — > 0) and again arrive 
at eq. Q determining no(p). The entropy S[n(p,T)] is given by the familiar expression [16] 

S[n(p,T)} = -2 |[n( P ,T)lnn( P ,r) + (l-n(p,T))ln(l-n(p,r))]^. (8) 

It follows from eq. (JBJ) that the entropy Snfl related to the special solution no(p) contains 
the temperature independent term So = »Snfl (T —> 0) ~ x(pf —Pi)/pF- Thus, the function 
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Fig. 1 — H — T phase diagram of CeCoIns . Right panel: Superconducting-normal phase boundary [10] 
is shown by the solid and dashed lines with the solid square showing the point where the supercon- 
ducting phase transition changes from the second to the first order. The dotted line is given by eq. 
I|10|l and represents the transition T*(H) from the Landau Fermi liquid behavior (LFL) with the T 2 
regime in p(T) to the T 2 / 3 one. The solid line given by eq. 11H represents the crossover T*{H) from 
the T 2//3 regime in p(T) [3] to the non-Fermi liquid (NFL) behavior with p(T) oc T. Experimental facts 
obtained from resistivity measurements are shown by the solid squares [3,4]. The left panel shows the 
magnetic field dependence of T 2 Landau Fermi liquid coefficients of charge A(H) oc {H — H c o)~ 4 ^ 3 
and heat B(H) oc (H — H c q)~ 4 ^ 3 transport with experimental data taken from ref. [3,4]. 



n (p) represents the special solutions of both BCS and LFL equations determining the NFL 
behavior of HF liquid with FC. Namely, contrary to conventional BCS case, the FC solutions 
are characterized by infinitesimal value of superconducting gap, A — > 0, while both k(p) and 
S remain finite and S = 0. In contrast to the standard solutions of the LFL theory, the special 
ones no(p) are characterized by the entropy S containing the temperature independent term 
So. At T — > both the normal state of the HF liquid with the finite entropy 5*0 and the 
BCS state with 5 = coexist being separated by the first order phase transition, where the 
entropy undergoes a finite jump SS = So- Due to the thermodynamic inequality, 5Q < T8S, 
the heat SQ of the transition is equal to zero making the other thermodynamic functions 
continuous. Thus, both at the FCQPT point and behind it there are no critical fluctuations 
accompanying second order phase transitions and suppressing the quasiparticles. As a result, 
the quasiparticles survive and define the thermodynamic properties of the HF liquid. 

On the basis of the above special solutions related to FC, we can explain the main universal 
properties of the H — T phase diagram of the HF metal CeCoIns shown in fig. The latter 
substance is a d - wave superconductor with T c — 2.3 K, while a field tuned QCP with a critical 
field of H c q = 5.1 T coincides with H C 2, the upper critical field where superconductivity 
vanishes [3-5]. We note that in some cases H c0 = 0. For example, CeRu2Si2 shows no 
magnetic ordering down to lowest temperatures [18]. Therefore, in our simple HF model H c q 
can be treated as a fitting parameter. Under the application of magnetic field H c q, CeCoIn5 
demonstrates the NFL behavior down to T = [11, 19]. It also follows from the above 
consideration that H c q ~ H C 2 is an accidental coincidence. Indeed, H C 2 is determined by g 
which in turn is given by the coupling of electrons with magnetic, phonon, etc excitations 
rather than by H c q. As a result, under application of a pressure which influences differently 
g and H c q, the above coincidence will be removed in agreement with facts [6]. At relatively 
high temperatures, the superconducting-normal phase transition in CeCoIns shown by the 
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solid line in the right panel of fig. ^ is of the second order [9, 10] and S and the other 
thermodynamic quantities are continuous at the transition temperature T C (H). Since H C 2 ~ 
H c q, upon the application of magnetic field, the HF metal transits to its NFL state down to 
lowest temperatures as seen from fig. ^ As long as the phase transition is of the second order, 
the entropy of the superconducting phase Ssc(T) coincides with the entropy SnflC? 1 ) of NFL 
state, 

S SC (T^T C (H))=S NFL (T^T C (H)). (9) 

Since Ssc(T — > 0) — > 0, eq. © cannot be satisfied at sufficiently low temperatures due to the 
presence of temperature- independent term So- Thus, in accordance with experimental results 
[9, 10], the second order phase transition converts to first order one below some temperature 
Tq(H). To estimate T (H), we use the scaling idea of Volovik (see ref. [20] for details), who 
derived interpolation formula for the entropy of a d - wave superconductor in a magnetic 
field H, while Snfl has been estimated in [21]. As a result, upon using eq. ©, we obtain 
T (H)/T C ~ 0.3. This point coincides pretty well with experimental value, shown on the fig. 
^ We note that the prediction that the superconducting phase transition may change its order 
had been made in the early 1960-s [22]. Being based on the thermodynamic considerations, 
our proof is robust and can be expanded on cases when the superconducting phase is replaced 
by another ordered state. Namely, if the superconducting phase were replaced by some other 
ordered phase separated from the NFL phase by the second order phase transition at H = 0, 
then at some temperature Tq(H) this phase transition should change its order. It follows 
from above consideration, that NFL phase has the temperature independent entropy term So- 
Since in the ordered phase the Nernst theorem (S — > as T — > 0) should hold, we conclude 
that there is the entropy step (from So to zero) as T — > while a system traverses the phase 
transition line from ordered phase to the NFL one. This means that this phase transition 
should change its order at Tq(H). For example, we predict that the AFM phase transition in 
YbRh 2 Si2 with Tn(H) (representing the field dependence of Neel temperature) should become 
first order at T < Tq(H), where Tq(H) is some finite temperature. Under constant entropy 
(adiabatic) conditions, there should be a temperature step as a magnetic field crosses the 
above phase boundary due to the above thermodynamic inequality. Indeed, the entropy jump 
would release the heat, but since S = const the heat is absorbed, causing the temperature 
to decrease in order to keep the constant entropy of the NFL state. Note that the minimal 
jump is given by the temperature-independent term So, which can be quite large so that the 
corresponding HF metal can be used as an effective cooler at low temperatures. 

The entropy Snfl determines the anomalous behavior of CeCoIn5 in the NFL region of 
the phase diagram. The term So ~ x{pf —Pi)/pF can be determined from the experimental 
data on spin susceptibility (following Curie law) and the specific heat jump AC at T c [21]. 
In HF metals like CeCoIns the normalized jump AC/C n ~ 4.5 is substantially higher than 
the ordinary BCS value [23], where C n is the specific heat of a normal state. The specific 
heat jump is not proportional to T c and is related to the fermion condensate parameter 
5pFC — {Pf—Pi)/PF ~ Sq/x, therefore the normalized jump AC/C n can be large [21,24]. This 
estimation gives SpFC — 0.044 [21]. The entropy Snfl determines also both thermal expansion 
coefficient a — —dS/dP and Griineisen ratio T = a/C n [21,25,26]. Since the entropy has the 
temperature independent part So, the thermal expansion coefficient a ~ —dSo/dP becomes 
temperature independent at low temperatures. Therefore, at T — > 0, a(T) — -> const, while the 
specific heat C n (T) — ^ 0. As a result, T(T — * 0) diverges in coincidence with the facts [11]. 

Now we consider the LFL behavior tuned by a magnetic field H > H c o- The LFL regime 
is characterized by the temperature dependence of the resistivity, p{T) = po + A(H)T 2 , with 
po being the temperature independent part and A(H) is the scattering coefficient. Since the 
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NFL behavior of CeCoIris coincides with that of YbRh2(Sio.95Geo,o5)2 and YbRh2Si2 [27,28] 
we would expect that the LFL behavior of these substances also coincide. For example, in 
YbRli2(Sio. 95660.05)2 the scattering coefficient diverges as A(H) cx (H — i/ c o) _1 [27] while in 
CeCoh-15 it diverges as A(H) cx (H — H c q) c with the exponent c ~ —4/3 [3-5]. In magnetic 
fields, the exponent c = — 1 characterizes the function A(H) of HF liquid with FC [19], while 
the exponent c = —4/3 describes the function A(H) of HF liquid on the disordered side 
of FCQPT [19,29,30]. To understand this striking change in the behavior of CeCoIns, we 
recall that FC has just appeared in this substance since Spfc — {Pf — Pi)/PF — 0.044 <C 1. 
As soon as the magnetic field is sufficiently high, H > H CT , (H CI is a critical field moving 
the HF liquid from the ordered side of FCQPT to the disordered side), Zeeman splitting 
6pp — (pfi~Pf2)/pf °f t ne two Fermi surfaces of HF liquid exceeds the condensate parameter, 
5pf > SpFC, an d the HF liquid with FC becomes LFL placed on the disordered side near QCP. 
Here ppi and pfi are the Fermi momenta of the two Fermi surfaces formed by the application 
of a magnetic field. The splitting can be estimated as p F SpF /M*(H) ~ HfiB, where \ib is the 
Bohr magneton. Taking into account that A(H) cx (M*(H)) 2 we obtain (H cr — H c q)/H c q ~ 
(ciSpf) 3 - Our estimations of the coefficient ci based on the experimental function A(H) 
show that ci ~ 5, and we obtain that the reduced field (H cr — H c o)/H c q ~ (ciSpp) 3 — 0.02. 
Thus, we can safely suggest that the reduced field of 0.02 is much smaller than the minimal 
reduced field 0.1 where A(H) measurements have been carried out in ref. [3,5]. As a result, the 
electronic system of CeCoIn 5 is placed on the disordered side of FCQPT by the application of 
such a high field and reveals A(H) cx {H — H c0 )~ 4 / 3 . We can see from the left panel of fig. 2] 
that the coefficient B(H) has the same critical field dependence. Here B(H) stands for the In- 
dependent contribution to the thermal resistivity and is related to A(H) by a field-independent 
factor, A{H) / B(H) ~ 0.47, as it should be in the case of ordinary metals [3,4] and HF metals 
demonstrating the LFL behavior. At sufficiently low temperatures and decreasing field when 
H < H cr , we predict that CeCoIns demonstrates the LFL behavior while the exponent c will 
change from c = —4/3 to c = —1. 

At low temperatures and H ~ H cri the system remains in the LFL regime, but at elevated 
temperatures there exists a temperature T*{H) where the influence of FC related to So is 
recovered and the NFL behavior is restored. To calculate the function T*(H), we note that 
the effective mass M* cannot be changed at T*(H). Since at T > T*(H) the effective mass 
M*(T) cx l/T [31] and at T < T*(H), M*(H) oc (H - # c0 )~ 2/3 , we have 



The function T*(H) given by eq. I|10|) is represented by the dotted line in the right panel of 
fig. H In high magnetic fields, H ^> H cr , there is a new crossover line because the effective 
mass starts to depend on temperature as M*(T) cx T~ 2 / 3 [29,30] and T*(H) becomes 



The crossover line given by eq. Ijll|l is represented by the solid line in fig. ^ (right panel) . As 
it is seen from fig. ^ the behavior of these lines described by eqs. IjlOII and (jl 1|> are in good 
agreement with experimental facts [3,4]. Eventually, when T/ > T > T*(H) the influence of 
FC determined by So restores and the system demonstrates the NFL behavior with M* cx l/T 
and p(T) oc T [19,30]. Here Tf is the temperature at which the influence of FC vanishes. For 
example, it can be estimated by using the condition, So <C S(Tf). The NFL behavior related 
to the So term can also be observed in measurements of tunneling conductivity and dynamic 
conductance which are expected to be noticeably asymmetrical with respect to the change of 
voltage bias from V to — V in HF liquid with FC [32]. Such asymmetrical conductivity was 



T*(H) cx {H-H c0 ) 2/3 . 



(10) 



T*(H)oc{H-H c0 ). 



(11) 
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recently observed experimentally in CeCoIns [33]. The behavior of the conductivity can be 
specific when the HF metal transits from its LFL state induced by the application of magnetic 
field to NFL one at elevated T. We predict that in the case of CeCoIns the conductivity being 
symmetrical in the LFL regime becomes gradually asymmetrical reaching its maximum in the 
NFL state at elevated temperatures when T > T*(H) and eventually vanishes. 

In summary we have presented for the first time theoretical description of the whole 
phase diagram of CeCoIns including the change of the second order superconducting phase 
transition to the first one under the application of rising magnetic field. We have shown that 
quasiparticles survive down to lowest temperatures. Our description of the HF metal CeCoIn 5 
based on the notion of quasiparticles and FCQPT is in good agreement with facts. 

We thank P. Coleman for stimulating discussions. This work was supported in part by 
RFBR, project No. 05-02-16085. The visit of VRS to Clark Atlanta University has been 
supported by NSF through a grant to CTSPS. 
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